Spin-orbit fields in ferromagnetic metal/semiconductor junctions 
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The existence of the spin-orbit fields induced by the interface structure in Fe/GaAs junctions is 
proven from first-principles calculations. While the underlying symmetry of the fields follows that 
of the interface, the specific realization of the symmetry depends on the electron momentum and 
energy. The calculated atomic-layer-resolved expectation values of the Bloch states' spins show that 
the spin-orbit fields peak at the GaAs side of the interface. The employed technique is applicable 
to ferromagnetic junctions in general. 
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To manipulate electron spins one usually resorts to 
the magnetic field. Recent experimental advances show 
that the relativistic spin-orbit interaction can do the job 
as well. Indeed, spin-orbit coupling appears as an ef- 
fective magnetic force acting on the electron spin and, 
equally important, on the electron motion These 
effective magnetic — more descriptively spin-orbit — fields 
emerge as a collective effect of individual atomic spin- 
orbit couplings in systems lacking space inversion sym- 
metry. By employing first-principles techniques, here we 
uncover the microscopic structure of the spin-orbit fields. 
We demonstrate the origin of their effective spectral de- 
scription widely used to study magnetic, transport, and 
optical manifestations of spin-orbit phenomena in solids. 
We visualize the charming patterns of the microscopic 
spin-orbit fields, atomic layer by layer, and show their 
interface origin for the technologically important case of 
Fe/GaAs nanojunctions recently used in tunneling mag- 
netoanisotropy experiments These spin-orbit fields 
give important microscopic clues in searching for verti- 
cal tunneling and lateral magnctotransport phenomena 
in single ferromagnetic layer metal/semiconductor junc- 
tions in which the interface, not the bulk, dictates the 
magnetic symmetry. 

In solid-state systems lacking space inversion symme- 
try spin-orbit coupling (SOC) is manifested on the spec- 
tral level as an effective magnetic field whose direction 
and magnitude depend on the electron momentum Q . In 
other words, spin-up and spin-down electrons of the same 
momentum have different energies. The most promi- 
nent examples are the Dresselhaus spin-orbit coupling 
[3| describing the effects of the bulk inversion asym- 
metry (BIA) as in zinc-blende semiconductors, and the 
Bychkov-Rashba spin-orbit coupling [||, describing the 
effects of the structure inversion asymmetry (SIA), as in 
asymmetric quantum wells. Apart from semiconductor 
structures, where the Bychkov-Rashba coupling has been 
extensively studied y, la Ml it has also been investigated 
in metallic surfaces [a El, I lOL 11, 12, EH, in graphene on 
a Ni substrate [li] ] , a nd in Au and Ag monolayers on 
W(110) substrates [H]. 



Here we go beyond the effective spectral description. 
We ask how the effective spin-orbit fields (SOFs) are 
distributed in real space, what is their typical magni- 
tude, functional form, and symmetry pattern. We ex- 
plore these issues on the example of Fe/GaAs junctions 
in which the phenomenon of tunneling anisotropic mag- 
netoresistance (TAMR) @, EB, Ell 

was observed [3 and 
explained in terms of the interface-induced spin-orbit 
fields. The magnetization direction in Fc makes for a 
nice control knob orienting the spin, while the mag- 
netic anisotropy is determined by the interface sym- 
metry, not by the symmetry of the bulk materials. 
These two features allow for investigating the micro- 
scopic nature of the spin-orbit fields. Such ferromagnetic 
metal/semiconductor junctions are useful for room tem- 
perature spin injection and magnetotransport HI, 13, 2oL 

In earlier studies of SOC effects in surfaces 



12l . Il3l . |23| and interfaces [FJ, ll5( the SOC parameters 
have been extracted by fitting the energy bands close 
to the r point, assuming a Bychkov-Rashba-typc SOC. 
This procedure, however, requires a priory knowledge of 
the specific form of the SOF and applies only to very 
small fc-vectors. Here we introduce a procedure to ob- 
tain the form of the SOF for any fc-point directly from 
ab-inito data. We illustrate this on a Fc/GaAs slab, in 
which the magnetization is perpendicular to the layers. 
This configuration allows for a simple relation between 
the spin-expectation values and SOF. By computing the 
spin expectation values from first principles, the specific 
vector pattern of the SOF, unknown a priory, is obtained. 
The results show a highly anisotropic SOF, which takes 
on different forms, from the ones known in semiconduc- 
tor physics to more exotic ones for states away from the 
r point. One fascinating outcome is the possibility of 
flipping the anisotropy axis of the SOF when going from 
one band to another, consistent with the bias-induced 
inversion of the TAMR observed in experiments 0, E3] • 
In order to deduce the analytical form of the SOC 
Hamiltonian at low electronic momenta in our Fe/GaAs 
slab, we propose a general model from which the SOF 
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can be extracted when contrasted with the ab-inito data. 
The non-centrosymmctric GaAs layer is of the Did sym- 
metry, thus exhibiting the BIA spin-orbit coupling. The 
interface lowers the symmetry to Ci v with the twofold ro- 
tation axis C2 along the growth direction [001] Q. The 
Civ symmetry accounts for both the BIA and SIA; the 
spin-orbit field lies in the plane of the slab, perpendicular 
to the growth direction. The most general SOC Hamil- 
tonian consistent with Ci v symmetry can be written for 
the in-plane momenta around the T point as 

^SO = &(k X , ky)k X (Jy + f3(k X , ky)ky<J X , (l) 

where k x and k y are the components of the in-plane wave 
vector k, a x and o~ y are the Pauli matrices, and x and 
y correspond to the diagonal [110] and [110] crystallo- 
graphic directions in GaAs, respectively. The functional 
parameters a and (3 are even in the momenta 

a(k x , k y ) = a<°> + a {1) kl + a {2) k 2 y + ... , 
^k x ,ky)=l3^+P^k 2 x +^kl + ... . (2) 

The values of the expansion parameters cS l \ (3^ (i = 
0, 1,2,.. .) determine the specific form of the SOF. For 
example if a (0) = -/? (0) (a (0) = /3 (0) ), H so reduces in 
the limit of small k = |fe| to the well known as Bychkov- 
Rashba [B| (linearized Dresselhaus Q) SOCs. By intro- 
ducing the SOF field 

(f3(k x , ky)ky \ 
a{k x ,k y )k x , (3) 
/ 

Eq.([T]) can be rewritten as Ti so — w(k) ■ cr, where cr is 
the vector of Pauli matrices. 

When the magnetization is perpendicular to the lay- 
ers the in-plane components of the spin are influenced by 
the SOC only. Since the exchange field dominates over 
the SOC, the spin is quantized largely along the mag- 
netization direction, and the expectation values of the 
transverse components of the spin (s) can be obtained 
by considering H so as a perturbation. First order per- 
turbation theory gives 

(s x )=w x /A ; (s v )=Wy/A, (4) 

where A is the exchange splitting energy and H = 1. 
Using these relations, we can determine the specific form 
of w, which provides an effective spatially averaged SOF 
of the system. 

The small lattice mismatch between Fe (2.87 A) and 
GaAs (5.65 A) allows for the smooth epitaxial growth 
of Fe on a GaAs (001) surface. Early investigations of 
the stability of 1 x 1 Fe/GaAs interfaces within density 
functional theory [25| revealed that when more than two 
atomic layers of Fe are deposited on a GaAs (001) sur- 
face, the flat or partially intermixed interfaces are more 



stable than the fully intermixed one. The large pd hy- 
bridization lowers the cohesive energy [26} and, there- 
fore, the antibonding bonds between Fe and As are more 
stable than the Fe-Ga bonds. While the density func- 
tional calculations indicate that the As-terminated flat 
interface is more stable than the partially intermixed one 
pjj ], a recent Z-contrast scanning transmission electron 
microscopy reported a single plane of alternating Fe and 
As atoms at an Fe/AlGaAs interface [27j • Here we con- 
sider an ideal structure with an As-terminated flat inter- 
face. Our slab, encased in a 6 A vacuum, contains 9 (001) 
atomic layers of GaAs with the diagonal lattice spacing 
d = a/\/2 = 3.997 A and three atomic planes of bec Fe. 

The electronic structure of the ideal Fe/GaAs slab 
has been calculated using the full potential linearized 
augmented plane waves (LAPW) technique [28| and a 
generalized gradient approximation for the exchange- 
correlation functional [29j. The SOC has been treated by 
solving the Dirac equation for the core electrons, while 
the valence electrons are treated within second varia- 
tional method. 

The band structure of the Fe/GaAs heterostructure 
along the high symmetry lines connecting the S — T — X 
points in the Brillouin zone (BZ) is shown in Fig. [IJ 
for a magnetization orientation along the [001] direction. 
Atomic-like contributions to the states coming from the 
interface atoms are emphasized by filled (red) and open 
(blue) circles for spin-up and spin-down, respectively. 




FIG. 1: Calculated band structure for the Fe/GaAs slab and 
magnetization along [100]. The states with spin-up (spin- 
down) at the Fe/GaAs interface are emphasized with red filled 
(blue open) circles whose radii are proportional to the corre- 
sponding charge density at the interface. The inset shows the 
As-terminated flat lxl interface model assumed in the study. 

We now explore the SOF in our heterostructure. Wc 
take bands (n = 1, . . . , 5) close to the Fermi level (see 
the band labels in Fig. []} and calculate the expectation 
values of the transverse spin components close to the T 
point. Thus by keeping only the linear terms in k of the 
SOF and fitting w x / A and w y /A to the ab-inito data for 
(s x ) and (s y ), respectively, we can determine the band- 
resolved Bychkov-Rashba-type, a = (a^ — /3^°- ) )/2, and 
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Dresselhaus-type, 7 = (a'- -' +/3^)/2, parameters. They 
arc shown in Tablc|TJ Interesting, the product changes 
sign from band n = 1 to n = 2, indicating a flipping of the 
symmetry axis of the SOF (compare left and right panels 
in Fig. [J). While typically the ratio a/7 is measured (30| . 
here by choosing A = 2 eV, which is within ±20% range 
of the exchange splitting of the Ti2 and T25' states in bec 
Fe, we obtain absolute values for a and 7. The extracted 
values are similar to those found in semiconductors [3J]. 



band label (n) 


1 


2 


3 


4 


5 


a/7 


-0.934 


1.161 


1.504 


0.805 


1.866 


a [eVA] 


-0.033 


0.125 


0.091 


-0.010 


0.090 


7 [eVA] 


0.035 


0.108 


0.060 


-0.013 


0.048 



TABLE I: Band-resolved Bychkov-Rashba-type and 
Dresselhaus-type spin-orbit coupling parameters. The 
states close to the Fermi level (see labels in Fig. [TJ and 
momenta close to the F point with k = 7r/100d are considered 
for fitting the values. The exchange splitting A = 2eV. 

We now concentrate on the two bands in vicinity of 
the Fermi level. The left panel in Fig. [5] corresponds to 
band n = 1, while the right one to the n = 2. The three 
different cases, k = n/50d, 7r/8d and 7r/5d, are considered 
for the plots of SOF, w(k) (bottom parts), and the polar 
plots of the field strength w = \w(k)\ (upper parts) in 
units of the exchange splitting A. Close to the T point the 
SOF resemble the interference of Bychkov-Rashba-type 
and Dresselhaus-type SOCs (see Fig. [2k). However, away 
from the V point higher in k terms become relevant and 
more exotic patterns in the SOF manifest (see Figs[2}3,c). 
Starting with k = 7r/50d, where linear terms in the SOF 
are dominant, the cubic terms in the SOF are already 
present for k = ir/8d, and they get more pronounced 
together with terms of fifth and higher orders for k = 
n/5d. Those terms we have to include for fitting the 
SOF field model (see curves in polar plots in Fig. [2]) to 
the ab-inito data (circles). The Ci v symmetry of the SOF 
is preserved for all k. 

In order to investigate from where the major contri- 
bution to the SOC in Fe/GaAs heterostructure comes, 
we calculate the layer-resolved SOF. In Fig. [3] we plot 
the Bychkov-Rashba-type and Dresselhaus-type parame- 
ters for each atomic layer within our slab for two bands, 
n = 1,2. The values have been obtained by fitting the 
linearized SOF to the ab-inito data. By studing the na- 
ture of the bands we find that the band n = 1 exhibits 
an interface character, while the band n = 2 is a surface 
state corresponding to the As-terminated (001) surface 
of the slab (see also band character in the vicinity of the 
r point in Fig. []}. In both cases the maximal values of 
the SOC parameters are near the Fe/GaAs interface, and 
especially, at the next nearest Ga atom plane as can be 
seen in Fig. [3] The maximal effect of the SOC is shifted 
to the GaAs region due to the fact that the interface As 
plane is in direct contact to the Fo and gets strongly po- 




FIG. 2: Spin-orbit-field "butterflies". Left panel corresponds 
to band n = 1 and the right one to n = 2. Polar plots of the 
spin-orbit couling fields as well as the corresponding vector 
fields are shown for a, k = 7r/50d; b, k = n/8d; c, k = 
ir/5d. The lines in a, and b, are fits of the spin-orbit field 
(in units of the exchange splitting) to the ab-initio data for 
|(s)|/fc shown by circles. The lines in c, are guides for the 
eye. The corresponding vector fields (s) are shown below the 
polar plots. The length of the vectors has been rescaled for 
convenience. 

larized due to the proximity effect. Consequently, the 
in-plane spin components, and therefore, the SOF com- 
ponents are considerably reduced. Figure [3J also reveals 
that the main contribution to the SOF comes from the 
interface even in the case of the surface band (n = 2). 
The SOFs may look different for different geometries, 
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FIG. 3: Atomic-layer resolved linearized Bychkov-Rashba- 
type and Dresselhaus-type spin-orbit fields. SIA and BIA 
parameters close to the F point, k — 7r/100d, are plotted as 
a function of the atomic layer in the slab for a, the interface 
band n = 1 and, b, the surface band n = 2. Parameters a 
and 7 are in units of eVA, using A = 2 eV; in the Fe layers 
they correspond to the sum of the spin moments of two Fe 
atoms within a single layer. 

growth directions and/or different materials composi- 
tions (say, Fe/GaAs/Au) but the methodology pre- 
sented here can be applied to rather generic cases. The 
SOF symmetry is crucial for fully understanding spin- 
injection from ferromagnets into semiconductors as well 
as magnetic and transport anisotropic properties in such 
junctions. This interface vector field can be observed 
by X-ray magnetic circular dichroism, while its manifes- 
tations are to be seen in magnetotransport experiments 
which effectively sample, with certain weights, the whole 
fc-space. 

The extracted magnitudes of a and 7 are within the 
range found in semiconductors 0]. However, the magni- 
tude of the SOC parameters can differ considerably from 
band to band. Similarly for the sign. This has observ- 
able consequences in transport properties of Fe/GaAs 
hetero junctions such as the TAMR Q. It has already 
been proposed that the Cn v symmetry of the TAMR in 
ferromagnet/scmiconductor-bascd tunnel junctions orig- 
inates from the interference between Bychkov-Rashba- 
type and Dresselhaus-type SOCs and that the sign of the 
TAMR is determined by the product a"/ (here a and 7 
refer to the average values of the corresponding SOC pa- 
rameters) @, [H . It has also been suggested that the sign 
of «7 (and therefore that of the TAMR) could be tuned 
by an external bias @,[2l|. This possibility is consistent 
with our results, in particular for the bands n = 1,2, 
where the product changes sign. 
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